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Abstract The SuperCDMS collaboration is a dark matter search effort aimed at 
detecting the scattering of WIMP dark matter from nuclei in cryogenic germa- 
nium targets. The CDMS Detector Monte Carlo (CDMS-DMC) is a simulation 
tool aimed at achieving a deeper understanding of the performance of the Super- 
CDMS detectors and aiding the dark matter search analysis. We present results 
from validation of the phonon physics described in the CDMS-DMC and outline 
work towards utilizing it in future WIMP search analyses. 
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1 Introduction 



The SuperCDMS phase of the Cryogenic Dark Matter Search (CDMS) experi- 
ment utilizes cylindrical germanium targets patterned with cryogenic Transition 
Edge Sensors (TESs) and ionization sensors in an effort to detect WIMP dark 
matter directly scattering from atomic nuclei 1 ' 2 . SuperCDMS detectors, known 
as iZIPs, consist of interleaved ionization sensors and TES phonon sensors with 
phonon-collecting aluminum fins patterned on both the top and bottom surfaces of 
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a 1" thick germanium crystal 3 ' 4 . The phonon signal is read by multiple channels 
to provide information about event position. The CDMS Detector Monte Carlo 
(CDMS-DMC) simulates phonon and charge carrier propagation as well as the 
TES and ionization sensors to provide a new tool for investigating the behavior 
of the CDMS detectors. Tuning and validation of the simulated phonon physics is 
necessary to ensure that the CDMS-DMC accurately reproduces detector behavior. 
This validation is performed by comparing results from CDMS-DMC simulations 
with data taken at CDMS testing facilities and at the underground experimental 
site. 



2 Distributions of Event Parameters 

In recent years, the CDMS collaboration has produced a number of new detectors 
with various TES sensor layouts on the surface. This variety of detector topologies 
provides an opportunity to validate the phonon physics in the CDMS-DMC across 
multiple detectors 5 - 6 , which prevents 'over-tuning' of physical parameters to ef- 
fects stemming from particular detector layouts and allows for confirmation of the 
validation procedure. Much validation work has focused on the CDMS-DMC's 
ability to reproduce distributions of event parameters (energy partitioning into the 
different TES channels, pulse rise and fall times, etc.) produced when the detec- 
tor is exposed to deeply-penetrating gammas from a calibration source. Figures 
[TJ and |2] show two example distributions from real data and simulations of events 
dispersed throughout the bulk of an iZIP detector with four TES channels on each 
surface, split into three inner and one outer channels. Figure [T] shows a histogram 
illustrating the partitioning of phonon energy into the 8 TES channels, while Fig- 
ure [2] shows a histogram of the phonon pulse decay times, which are determined 
by how rapidly phonons are removed from the crystals by the phonon-collecting 
Al fins. Comparing these distributions and many others has allowed us to tune the 
probability of phonon absorption vs. reflection at the patterned surfaces, compare 
measurements of the TES critical temperatures, consider different values for the 
anharmonic decay constant (eventually converging to near the literature value), 
and validate the area of aluminum coverage on the detector surface used in the 
CDMS-DMC. 



3 Phonon Pulse Matching 

The next phase of validation of the CDMS-DMC involves directly matching the 
sets of TES pulses observed in real data events with sets of pulses from events 
simulated in the Detector Monte Carlo. Matching TES pulses observed in data 
to simulated pulses requires not just accurate phonon physics but also a well- 
tuned description of the TES parameters, which is described elsewhere in these 
proceedings 7 ' 8 . The procedure for matching the TES pulses begins with thousands 
of electron recoils simulated on a grid of points throughout the germanium crystal. 
A single data event is selected, and the set of phonon pulses is compared to the 
phonon pulses from each of the simulated events. Each simulated event is fit to 
the data event using a single amplitude scaling parameter. A weighted sum of the 
residuals between the data pulses and the simulated pulses is used as the goodness 
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Fig. 1 (Color online) Distribution of energy 
partitioning (energy in one channel / total event 
energy) in real data events (solid blue) and 
CDMS-DMC simulations (dashed red) of an 
iZIP4 detector. The detector has 8 total chan- 
nels, 4 on each detector surface, divided into 
3 inner channels and 1 outer channels. Each 
channel covers 1/4 of the detector surface. Each 
event contributes 8 entries to the histogram, 
one for each channel. The y-axis represents the 
counts per bin normalized to sum to 1. The 
difference in widths indicates that the data ex- 
hibits slightly stronger energy partitioning than 
the simulation. 
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Fig. 2 (Color online) Decay times of TES 
pulses observed in data (solid blue) and in 
CDMS-DMC simulations (dashed red). The de- 
cay of each TES pulse is fit to a single falling 
exponential from 750 to 2750 [is after the onset 
of the event. The decay time distributions are 
then fit to Gaussians to provide the means and 
widths indicated in the legend; these fits are not 
shown. The mean decay times of data and sim 
are in agreement; the differences in the widths 
of the distributions are most likely due to chan- 
nel variations and noise. The secondary peak in 
the simulation histogram is due to a single TES 
channel. The y axes used for the two histograms 
are scaled relative to one another to allow for 
comparison. 



of fit metric; after a data event has been compared to all simulated events, the 
simulated event with the best (lowest) goodness of fit metric is selected as the 
best match for that particular data event. Figures [3] and |4] show the results of this 
matching routine for a single illustrative data event and shows that the selected 
simulated event matches the data pulses quite closely. Figure [3] shows the entire 
sampling length of the TES pulses, w2850 jJ.s, while Figure [4] zooms into the first 
400 jis of the pulse, which contains most of the position-dependent variation in 
pulse shape. 

The ultimate goal of the TES pulse matching work is to demonstrate pulse- 
shape based discrimination of background electron recoils in the bulk and surface 
regions of the detector from the nuclear recoils that represent the WLMP-nucleon 
scattering signal. Figure [5] shows the weighted sum of residuals when bulk elec- 
tron recoils in the data are compared with simulated bulk electron recoils and sim- 
ulated bulk nuclear recoils. At present, these distributions are not well-separated, 
and in fact many of the electron recoils from the data find better matches among 
the simulated nuclear recoils than the simulated electron recoils. Thus, this pulse 
matching analysis does not yet demonstrate the ability to discriminate signal from 
background based solely on simulated pulse shapes. 
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Fig. 3 (Color online) Set of phonon pulses from the 8 TES channels of the iZIP4 detector from 
a single data event (solid blue) shown with the best match found among the set of simulated 
events (dashed red). The data was taken at the UC Berkeley SuperCDMS testing facility. The 
left and right columns show pulses from the 4 channels on the top and bottom surfaces of the 
detector, respectively. Channel A is the outer TES channel, and channels B, C, and D are the 3 
inner channels on each side. The y-axis represents power in the TES in arbitrary digitized units. 
The black vertical lines indicate the range (400 fis) of the zoomed figure shown below. 




Fig. 4 (Color online) Zoomed version of figure above, showing a set of phonon pulses from a 
single data event (solid blue) shown with the best match found among the set of simulated events 
(dashed red). Pulses are zoomed in to the first 400 /is. 



5 




Fig. 5 (Color online) Distribution of the sum of weighted residuals (quality of pulse match) for 
bulk electron recoils from the data matched with simulated bulk electron recoils (solid blue) and 
simulated bulk nuclear recoils (dashed black). There is a complete lack of separation between the 
two distributions, and interestingly, the bulk NR simulations seem to provide better matches to 
real electron recoils than the simulated electron recoils do. The y-axis is the counts per histogram 
bin normalized to sum to 1. 



4 Future Work 

Validation of the CDMS-DMC is an ongoing process of continuous refinement. 
Additional work regarding phonon transport Monte Carlo simulation can be found 
in the proceedings 9 . Future CDMS-DMC work aims to create a model-based 
WIMP search analysis method utilizing the pulse matching scheme presented 
above. A first step toward this goal involves quantifying the degree of observed 
variation in the phonon pulse shapes as the event position changes. These results 
can then be used to determine how reliably simulated pulses at different posi- 
tions can be interpolated to achieve stronger agreement between simulation and 
data. The simulations of electron and nuclear recoils must also be independently 
compared to the data to determine whether the CDMS-DMC produces the excess 
'peakiness' of electron recoils demonstrated in other work 10 , as this feature will 
be critical to a WIMP search analysis based on pulse-shape discrimination using 
the CDMS-DMC. 



5 Conclusions 

The CDMS-DMC is able to reproduce many features of the phonon and TES be- 
havior observed in real CDMS detectors. This agreement has been demonstrated 
using broad distributions of multiple event features such as energy partitioning and 
pulse shapes and also on event-by-event matching of sets of real TES pulses from 
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CDMS detectors to pulses simulated using the CDMS-DMC. The goals of match- 
ing individual events from real SuperCDMS data with events simulated using the 
CDMS-DMC is to understand the position-dependent variation in the pulse shapes 
and to perform a WIMP search analysis using pulse-shape-based discrimination 
of background electron recoils from the nuclear recoil signal. 
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